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ABSTRACT 

We present results from dynamical Monte Carlo simulations of dense star 
clusters near the Galactic center. These clusters sink toward the center of the 
Galaxy by dynamical friction. During their inspiral, they may undergo core 
collapse and form an intermediate-mass black hole through runaway collisions. 
Such a cluster can then reach within a parsec of the Galactic center before it 
completely disrupts, releasing many young stars in this region. This scenario 
provides a natural explanation for the presence of the young stars observed near 
the Galactic center. Here we determine the initial conditions for this scenario 
to work and we derive the mass distribution of cluster stars as a function of 
distance from the Galactic center. For reasonable initial conditions, we find 
that clusters massive enough for rapid inspiral would include more massive stars 
(m^ > 30 Mq) than currently observed in the inspiral region. We point out 
several possible explanations for this apparent discrepancy. 

Subject headings: black hole physics — Galaxy: center — galaxies: star clusters 
— methods: N-Body simulations, stellar dynamics 



1. Introduction 

Multi- wavelength observations of the Galactic center (GC) have revealed a supermassive 
black hole (SBH) and a stellar cusp surrounding it. The mass of the SBH, Msbh — 4x10^ Mq, 
is one of the most reliable estimates for massive black holes at the centers of galaxies. The 
stellar population within Ipc of the SBH contains a variety of young and massive stars. 
Some of them are only about 20 Myr old and get as close as a few light days to the SBH; 
while from 0.1 to 0.4 pc even younger stars are found with ages 3 to 7 Myr. The presence of 
these young stars in the immediate vicinity of the SBH poses a problem known as the youth 
paradox. Their in situ formation is problematic since the SBH has a strong tidal influence 
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in this region. However, the time required for the migration of these stars from > 1 pc by 
dynamical friction would exceed their inferred ages, unless the migration rate is somehow 
accelerated. 

In addition to their youth, these stars exhibit some peculiar dynamical properties. One 
of them is strong clustering. The sources IRS 16 and IRS 13 were initially thought to be 
single sources but later resolved into multiple components. IRS 16 is a collection of young 
(3-7 Myr) He I emission-line stars lying l"-7" in projection from Sgr A*. They form a co- 
moving group that is counter-rotating with respect to the Galaxy. IRS 13 is another complex, 
composed of hot bright stars, about 3.6" south-west of Sgr A*. In particular, IRS 13E is a 
compact, massive star cluster which is a few million years old (Maillard et al. 2004). Genzel 
et al. (2003) and Horrobin et al. (2004) find that the young He I emission-line stars in 
the vicinity of Sgr A* are all concentrated in two disks, forming two co-moving (but not 
gravitationally bound) populations. In addition to the clustering, stars within ~ 0.04 pc of 
Sgr A*, which are about 10 Myr old, seem to have higher than normal eccentricities (Schodel 
et al. 2003), but this anomaly may be explained by selection effects (Ghez et al. 2004). 

One possible explanation for the presence of these young stars is that they are not young 
but rejuvenated. A possible path to rejuvenation, stellar mergers, is shown not to be viable 
by Genzel et al. (2003). Another path, squeezars, heating of stars by close tidal encounters 
with the SBH also fails to explain the observed population (Alexander & Morris 2003). 

Aside from rejuvenation, another possible scenario, especially suggested by the obser- 
vations of the two counter rotating and coeval disks of young stars around the SBH, is the 
infall and collision of two molecular clouds in the vicinity of the SBH. Such a collision may 
provide the required densities for star formation. Genzel et al. (2003) suggest that the stars 
resulting from the two colliding clouds will form two co-moving populations. However, since 
a collision leading to high densities will be highly inelastic, it is not clear that two distinct 
populations will form. To our knowledge, the collision and further evolution of dense clouds 
in the vicinity of a SBH has never been studied in detail. So, it is hard to decide whether 
this is a viable scenario. 

A possibility that we investigate in this paper is the inspiral of a star cluster by dynamical 
friction, as suggested by Gerhard (2001). Portegies Zwart et al. (2003) carried out iV-body 
simulations to study this scenario. The clusters in their simulations disrupted at > 1 pc 
from the GC, so cannot explain the presence the young star population very close to Sgr A*. 
Kim & Morris (2003), again using A^-body simulations, concluded that, for a cluster to reach 
within the central parsec, it must either be very massive (> 10^ ^o) or have formed near the 
GC (at < 5pc). They found that, in both cases, a very high central density (~ 10® Mopc"^) 
is required, and concluded that this scenario is implausible. 



3 



Hansen & Milosavljevic (2003) suggested that the presence of an intermediate-mass 
black hole (IMBH) can stabilize a cluster so that survival all the way to within the central 
parsec can be achieved with much lower central densities. In addition, they suggested that 
an IMBH-SBH binary can perturb the young stars into radial orbits with small semi-major 
axes in a way similar to the Jupiter-Sun system creating short period comets (Quinn et al. 
1990). 

Kim et al. (2004) performed A^-body simulations similar to those of Kim & Morris (2003) 
but with an additional IMBH embedded in the inspiraling cluster. They found that an IMBH 
does decrease the requirement for a high central density, but its mass must be about 10% of 
the total cluster mass. Since this is much larger than estimates of the collapsed core mass in 
dynamical simulations (0.1-0.2%; see Portegies Zwart & McMillan 2002; Giirkan et al. 2004) 
they concluded that a realistic IMBH cannot help transport young stars into the central 
parsec. However, after core collapse, the central object can continue to grow by colliding 
with the stars that migrate to the center by relaxation. This growth continues until the 
massive stars start evolving off the main sequence and lead to cluster expansion by mass loss 
through winds and supernovae (Portegies Zwart et al. 2004; Freitag et al. 2004). 

Using our dynamical Monte Carlo code, we have carried out simulations where we form 
an IMBH through runaway collisions following core collapse during the cluster's inspiral 
towards the GC. We discussed some initial results in a previous paper (Giirkan & Rasio 
2004). Here we present all our results from simulations of dense star clusters near the 
GC using a more realistic Galactic mass distribution. We investigate the initial conditions 
required for these clusters to reach the GC within 3-10 Myr and undergo core collapse 
before disruption. In §2 we describe our numerical technique for cluster simulations and the 
inspiral mechanism. We present some semi-analytic calculations in §3, and the results of our 
full simulations in §4. We discuss these results and present our conclusions in §5. 

2. Initial Conditions and Numerical Technique 

For our simulations, we use a Monte Carlo technique, which provides an ideal compro- 
mise in terms of speed and accuracy. It has the star-by-star nature of the A^-body techniques 
but incorporates physical assumptions which lead to simplifications, reducing the computa- 
tion time. As a result, we can carry out simulations of systems with A^ > 10^ stars in < 100 
CPU hours. For details of our method, we refer to our previous work (see Joshi et al. 2000, 
2001 and references therein for basics, and Giirkan et al. 2004 for treatment of the realistic 
mass functions). Here we only explain in detail the additions to our code for this paper. 
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2.1. Initial conditions and Units 

For the initial mass function (IMF) we implemented a Salpeter distribution, dN/dm cx 
with mmin = 0.2 and mmax = 120 M0. Other choices of IMF, e.g., Miller-Scalo 
or Kroupa, do not change the core collapse time, as long as rrimax/'^min ^ 100 (Giirkan et al. 
2004). 

For the initial structure of our clusters we chose King models (see Binney & Tremaine 
1987, Chapter 4). This choice allows us to implement the initial tidal cutoff in a natural way. 
In addition, the rate of evolution of the system, which is determined by the central relaxation 
time (Giirkan et al. 2004), can be adjusted by a single parameter, the dimensionless central 
potential Wq. For the sake of simplicity we chose the tidal radius of the King model equal 
to the Jacobi radius of the cluster, (see §2.3) 

Throughout this paper we use standard Fokker-Planck units (Henon 1971): G = Mq = 
—4:Eq = 1, where G is Newton's gravitational constant, Mq is the initial total mass of 
the cluster, and Eq is the initial total energy. The conversion to physical units is done by 
calculating the physical mass of the system, and identifying the tidal radius of the King 
model used with the Jacobi radius. Various length and time scales of King models are given 
in Table 1 of Giirkan et al. (2004). 



2.2. Boundary Condition at the Center 

In the continuum limit, when a cluster undergoes core collapse, the central density be- 
comes infinite in a finite time. In addition to being unphysical, high densities require small 
timesteps and render the dynamical evolution of the system impossible to follow numeri- 
cally. When various physical processes and the finite radii of stars are taken into account, 
the central density cannot increase indefinitely. This can result from energy generation by 
"three-body binaries" (Giersz 2001, and references therein), "burning" of primordial binaries 
(Fregeau et al. 2003, and references therein) or physical collisions. For the young dense clus- 
ters we consider here, we expect energy generation from binary formation or binary burning 
to play a minor role, since for these systems, an interaction with a hard binary is likely to 
lead to a merger (Fregeau et al. 2004). 

The local collision time, i.e., the average time after which a star has experienced one 
collision, is given by 

^ n 30kms-^ R. M, ' ^ ' 
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where is the velocity dispersion, n is the number density of stars, fZ* and M* are the radius 
and the mass of the star under consideration. In the systems we consider in this paper, tcou 
is larger than the lifetime of the cluster except for the stars participating in the core collapse. 
We used an approximation to exploit this property. Rather than treating collisions explicitly 
as in Freitag & Benz (2002), we introduced a simple boundary condition near the center of 
the cluster: when a star is part of the collapsing core, we add its mass to a growing central 
point mass and remove it from the simulation. This central point mass is then used only 
during the calculation of the cluster gravitational potential. We determine whether a star is 
part of the collapsing core by monitoring its orbit's apocenter. We have chosen this criterion 
rather than one based on pericenter or instantaneous position to guarantee that the stars 
under consideration are restricted to a small region near the center. The threshold value we 
have chosen for the apocenter distance is 2 x 10~^ in Fokker-Planck units. For our models, 
tcow within this radius is ~ 10^ yr, i.e., all stars restricted to this region will rapidly undergo 
collisions. We determined the threshold value empirically. Choosing a large value leads 
to removal of stars from the system before core collapse. Choosing a very small value is 
also problematic because of the very short collision times. If collisions are not taken into 
account, a few massive stars sinking to the center of the cluster could easily reach energies 
comparable to the total energy of the cluster and have a substantial effect on the evolution 
of the system (Henon 1975, and references therein). We found that threshold values in the 
range 5 x 10~^-3 x 10~^ generally avoid these problems and 2 x 10~^ is a number suitable 
for all the clusters we simulated. 

We compared the rate of growth of the central mass in our simulations with results 
from other Monte Carlo simulations where collisions between the stars are treated more 
realistically Freitag et al. (2005,b). We found that we typically underestimate the rate of 
growth slightly, by 20 to 30%. We have also compared the evolution of the average mass 
among the innermost 5 000 stars and found that, in our simulations, the increase in average 
mass slightly lags behind, but reaches the same final value. This is probably related to the 
slower growth of the central mass in our code. It is not possible to determine which results 
are more accurate as these simulations cannot be repeated with more direct methods. In 
any case we do not expect this small difference to affect our results significantly. 



2.3. Tidal Truncation 

During the course of our simulations, we remove the stars that gain positive energies 
because of interactions, or whose apocenter lies beyond the cluster's tidal radius. This tidal 
radius depends on the mass distribution in the GC region and the cluster's current position. 
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The tidal radius of the cluster can be estimated by using the following expression for the 
Jacobi radius 

nij 



rj 



1/3 

R, (2) 



2M{R)_ 

where mj is the bound mass of the cluster, R is the distance from the GC, and M{R) is 
the mass within a circular orbit at this radius. The latest estimate for the Galactic mass 
distribution near the central BH is given by Genzel et al. (2003). They estimate the stellar 
mass density broken power law, 

p^{R) = 1.2 X 10' (^^^ MQpc-^ (3) 

with a = 1.63 for R ^ Rb, a = 1.0 for R > Rb, and Rb = 0.38 pc. This mass distribution is 
a self-consistent description that goes all the way from 0.1 to about 10 pc, but much farther 
than that it may underestimate the stellar density (Tal Alexander, private communication). 
We derive the formulae for dynamical friction on a point mass, resulting from a broken power 
law mass density distribution and an additional central BH of mass Msbh = 4 x 10' Mq 
(Melia & Falcke 2001; Ghez et al. 2003; Schodel et al. 2003), in Appendix A. As a result 
of the dynamical friction, the distance of the cluster from the GC changes continuously; we 
include the effect of this change by adjusting rj appropriately in our simulations. 

In principle, the treatment we use for dynamical friction is valid only for point masses 
and extended clusters require a modification of this treatment. However, the finite size 
affects only the Coulomb logarithm in equation A5, and, as long as the size of the cluster 
[vh) is small with respect to the size of the region that contributes to the dynamical friction, 
this leads only to a small decrease in the drag on the cluster (Binney & Tremaine 1987, 
§ 7.1). Since this condition is always satisfied in our simulations, we do not expect our 
results to be affected by this approximation. McMillan & Portegies Zwart (2003) established 
numerically the validity of the formulation we use by making comparisons to more exact 
A^-body simulations. 



3. Semi-analytical calculations 

For a point mass, the time required to reach the GC from an initial distance Rq can be 
calculated by integrating equations (A9) and (All). Since equation (A9) has a logarithmic 
singularity near the origin, for the calculations below, we stop the integration at i? = 0.01 pc. 
Most of the mass loss from the cluster takes place near the center, so neglecting the mass 
loss provides a reasonable first estimate for the inspiral time, t-m- Possible upper limits for 
tin are 3 or 10 Myr, which are about the lifetimes of the brightest IR stars observed near the 
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GC. By requiring the cluster to undergo core collapse before reaching the GC, we can also 
obtain a lower limit for ti^. This can be transformed into a condition on the initial structure 
of the cluster using t^c — 0.15ti.c(0) (Giirkan et al. 2004), where trc(O) is the initial central 
relaxation time. We illustrate these constraints in Figure 1. The solid lines correspond to the 
upper limits on tin from the stellar evolution timescale. Above these lines it takes more than 
the denoted time (3 or 10 Myr) to reach 0.01 pc for a point mass. The bow-like dashed lines 
correspond to lower limits for various King models. To the right of these lines, the inspiral 
time for a point mass is shorter than the core collapse time, hence runaway collisions cannot 
happen unless the system is initially coUisional. For these calculations we have computed the 
timescales for the King models by assuming their tidal radius is initially equal to their Jacobi 
radius given in Equation 2. In this figure we also show results of some of our simulations 
(see next section). 

The next step in approximation is using a realistic structure for the cluster but neglecting 
the changes in structure during the inspiral. Since the structure changes on a relaxation 
timescale and in the outer parts of the cluster the relaxation time is much longer than the 
inspiral time, this approximation provides a good estimate for the mass loss during the 
inspiral, at least until close to disruption. We have computed the rate of inspiral of different 
King models for different initial distances from the GC and initial cluster masses. We present 
our results for initial distance i?o = 10 pc in Figure 2. During the inspiral we calculate the 
change in bound mass as follows. First we calculate the tidal radius by using the current 
mass and distance from the GC by equation (2). From this and the structure of the given 
King model a new bound mass is calculated. This procedure is repeated until the bound 
mass converges to a value. Failure of convergence marks the disruption. 

In Figures 2 and 3 we show the results of our calculations for extended clusters. Unlike a 
point-mass, an extended cluster disrupts at a finite distance from the GC, -Rdis- This distance 
decreases as more of the cluster's mass is concentrated near its center, but is independent of 
the total mass in the cluster. The inspiral time, on the other hand, has a strong dependence 
on the total mass. The independence of i?dis from total mass can be understood as follows. 
For a given Wq, the disruption happens when the tidal radius reaches a specific fraction 
of the initial Jacobi radius. The relation between this value and the distance from the 
GC is independent of the initial cluster mass, since more massive clusters start with a 
proportionately larger initial Jacobi radius. This is also the main reason for the dependence 
of -Rdis on initial distance. Note that, in principle, a cluster can underfill its Jacobi radius 
at formation, i.e., clusters formed at large distances can have the same size as the clusters 
formed close to the GC. The disruption radii we obtain in this approximation are upper 
limits since mass segregation, the formation of a cusp, and the subsequent IMBH will all 
increase the mass concentration towards the center. 
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Fig. 1. — Region of Ro-ttiq space that allows fast inspiral (denoted by solid lines) and core 
collapse happening during inspiral (for various King models; denoted by dashed lines), in 
the point-mass approximation. Results of various simulations are shown with symbols of 
different shapes, corresponding to different initial structure; full symbols are for systems that 
experienced core collapse and open symbols correspond to systems that disrupted before core 
collapse. All models shown have tdis < 3 Myr except for the ones with gray symbols (colored 
red in the electronic edition), which have 3 Myr < t^is < 4 Myr 
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Fig. 2. — Rate of inspiral for point-mass and King models. The two families of curves 
correspond to initial masses mo = 3 x IO^Mq and mo = lO^M©. For each initial mass, 
results are shown for point-mass (solid line) and King models with Wq = 9 (long-dashed 
line), Wq = 6 (dashed line), and Wq = 3 (dotted line). 
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Fig. 3. — Dependence of disruption radius -Rdis on initial concentration and initial distance. 
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These calculations show that for a wide region of Ro~mo parameter space it is plausible 
that a cluster can reach to within a fraction of a parsec of the GC in a few million years. 
A high initial concentration [Wq > 8) is required for the cluster not to disrupt too far from 
the GC, as well as to guarantee rapid core collapse. More detailed simulations are required 
to verify and establish these findings, and to obtain the demographics of the stars that 
reach the GC. To illustrate this point, we present a comparison of point-mass and extended 
cluster models with a full simulation (Model5c, see next section) in Figure 4. All calculations 
presented in this figure start with the same initial mass and initial distance from the GC. 
This comparison shows that the effects of dynamical evolution, namely the enhanced mass 
loss due to expansion, can be quite important. Note that, for clusters with lower initial 
concentrations, the relaxation will be slower and the difference between full simulation and 
the extended static model will be smaller. However the deviations from the point-mass 
approximation will be larger. 



4. Results of simulations 

4.1. Overview of results 

We have performed simulations starting with King models of dimensionless central po- 
tentials Wq = 6, 8, and 9 at a variety of distances and initial masses. We present our results 
in Table 1. We ended our simulations when the number of stars in the cluster dropped to 
0.5% of its initial value. We defined this as the point of disruption and denoted the time to 
reach this point and distance from GC with t^is and Rdis, respectively. Allowing the cluster 
to evolve beyond this point will not change our results since, at this point, the mass is so 
low that dynamical friction is no longer effective, and the distance from the GC remains 
constant. In addition, as a result of the low number of stars, relaxation and the consequent 
evaporation of the cluster will be very rapid. Most of the clusters with high initial concen- 
tration went into core collapse before disruption and the central point mass described in 
§2 grew during the inspiral. This accumulated mass is indicated in column 7 of our table. 
For less dense clusters, the inspiral is faster than the core collapse and in those cases no 
central point mass is grown. We did not observe significant post-core-collapse expansion in 
our simulations. The growth of the central mass indicates that segregation of massive stars 
continue beyond the core collapse. 

An overview of the results shows that neglecting the changes in the structure of the 
cluster during its inspiral provides a good estimate for the time to disruption, t^is- However, 
this is not true for -Rdis? which is dependent on the total initial mass for an evolving cluster. 
This dependence is a result of the interplay between inspiral time and relaxation time. For 
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Fig. 4. — Comparison of a full simulation with an extended static model and the point-mass 
approximation. 
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Table 1. Results of Simulations 



model 


mo 


Ro 


Wn 








ID 


{Mq) 


(pc) 


(pc) 


(Myr) 


{Mq) 


la 


3 X 10^ 


20 


8 


1.04 


1.04 




2a 


3 X 10*' 


10 


8 


0.35 


0.31 




3a 


3 X 10^ 


5 


8 


0.22 


0.12 




4a 


z X iU 


on 


Q 

o 


l.Oo 


i.oy 


oUUU 


5a 


2 X 10^ 


10 


8 


0.88 


0.46 


580 


6a 


2 X 10^ 


5 


8 


0.34 


0.17 




7a 


10^ 


20 


8 


1.89 


3.16 


3600 


8a 


10^ 


10 


8 


0.98 


0.94 


2500 


9a 


10^ 


5 


8 


0.47 


0.32 


2000 


4b ... . 


2 X 10^ 


20 


6 


1.44 


3.00 




5b 


2 X 10^ 


10 


6 


0.42 


0.85 




6b 


2 X 10^ 


5 


6 


0.24 


0.31 




7b ... . 


10® 


20 


6 


2.09 


3.08 




8b ... . 


lO*" 


10 


6 


0.92 


0.88 




9b ... . 


10® 


5 


6 


0.42 


0.30 




10b . . . 


5 X 10^ 


20 


6 


4.59 


2.96 


560 


lib . . . 


5 X 10^ 


10 


6 


1.55 


0.86 




12b . . . 


3 X 10^ 


20 


6 


4.44 


3.00 


1200 


13b . . . 


3 X 10^ 


10 


6 


2.31 


0.80 




Ic 


7^ 

3 X 10® 


20 


9 


1.23 


1.10 


5500 


2c 


3 X 10® 


10 


9 


0.56 


0.35 


4500 


3c 


3 X 10® 


5 


9 


0.28 


0.14 


3000 


4c 


2 X 10® 


20 


9 


1.20 


1.70 


6000 


5c 


2 X 10® 


10 


9 


0.60 


0.52 


4900 


6c 


2 X 10® 


5 


9 


0.33 


0.20 


4200 


7c 


10® 


20 


9 


1.89 


3.23 


4700 


8c 


10® 


10 


9 


0.76 


1.04 


3500 


9c 


10® 


5 


9 


0.43 


0.34 


3600 


10c. .. . 


5 X 10^5 


20 


9 


3.14 


6.50 


2900 


11c. . . . 


5 X 10 '^ 


10 


9 


1.26 


1.95 


2700 


12c. .. . 


3 X 10= 


20 


9 


4.18 


11.0 


1900 


13c. .. . 


3 X 10= 


10 


9 


1.91 


3.04 


2000 


17c. . . . 


2 X 10® 


30 


9 


1.96 


3.54 


6800 


18c. .. . 


2 X 10® 


60 


9 


5.17 


13.8 


6800 


19c. .. . 


2 X 10® 


50 


9 


3.81 


9.63 


6600 



Note. — The first column is the model ID. The following three 
columns indicate the initial conditions of the cluster: mass in solar 
masses, initial distance from the GC in parsecs and the dimension- 
less central potential. The last three columns are the results of the 
simulations: the disruption distance in parsecs, the time it takes 
for disruption in Myr, and the accumulated mass at the center of 
the cluster in solar masses. If the cluster disrupts before going into 
core collapse, no mass is accumulated at the center. In these cases 
a hyphen is used in the last column. 
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a massive cluster containing a large number of stars, the inspiral time is short but the 
relaxation time is long. Hence only the central part of the cluster has time to evolve. For a 
less massive cluster the outer parts have time to evolve and expand. As a result, the mass 
loss is enhanced and the disruption happens at a further distance from the GC. 

In Figure 1 we plot the results of some our simulations to compare with point-mass 
estimates. We have only plotted models where disruption distance from the GC was small. 
For models started with Rq < 10 pc, only the ones within 1 pc of the GC; and for the others 
only the ones disrupted within 2pc of GC are plotted. We have used black symbols for 
models with tdis < 3 Myr, and gray symbols (colored red in the online version) otherwise. 
We find that inspiral and/or disruption time in simulations exceed the estimated inspiral 
time from point-mass approximation. Similarly, the disruption distance is always larger 
in full simulations than static cluster approximations. Finally, sometimes clusters disrupt 
without going into core collapse in simulations where the point-mass approximation predicts 
core collapse times shorter than inspiral times. These results show that semi-analytical 
approximations provide necessary but not sufficient conditions for our scenario to work. 
Increasing the realism by taking the tidal mass loss of the cluster and its expansion into 
account cannot be compensated by the increase in central density of the cluster, and leads 
to disruption happening earlier and/or further from the GC than indicated by the semi- 
analytical methods. It is not possible to resolve the core-collapse criterion, indicated by 
dashed lines in Figure 1, for Wq = 6, 9 models as is done for Wq = 8 case. For Wq = 6, the 
boundary is at low initial mass which correspond to low number of stars. Such clusters have 
short relaxation times and hence expand and disrupt rapidly before going into core collapse 
unless they start far from the GC. In contrast, the boundary for Wq = 9 lies at high initial 
mass and large number of particles, making the simulations of such clusters computationally 
impractical. 

Another trend that can be seen in Table 1 is that more mass is accumulated at the 
center for clusters that start at larger distances from the GC. This is simply because there 
is more time available for the mass buildup. This does not necessarily mean that a heavier 
IMBH will be formed in these systems, as the fate of a very massive star being bombarded 
by other stars is not certain. Also, being the total available mass, this is strictly an upper 
limit to rather than an estimate of the IMBH mass. 

In the rest of this section we present an analysis of the evolution of our Model 5c 
(see Figure 5), which started as a Wq = 9 King model, with mass mg = 2 x 10^ Mq at 
distance Rq = 10 pc. This cluster went into core collapse very rapidly (> 10^ years) and 
reached the GC in about 5 x 10^ years. After core collapse the central mass grew to a value 
Mcen — 5OOOM0, which is consistent with the estimate, Mcen — 0.002Mtotai, by Giirkan et al. 
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(2004). The disruption of the cluster took place at -Rdis = 0.6 pc from the GC. 

4.2. Stars Stripped from the Cluster 

During the inspiral, stars become unbound from the cluster as a result of the shrinking 
Jacobi radius. We present the mass distribution of these stars with respect to their distance 
from the GC in Figure 6, and the number of massive stars that leave the cluster within 
each distance bin in Figure 7. In these figures it is seen that the massive stars leave the 
cluster predominantly near disruption. This is because of mass segregation, but also because 
most stars leave the cluster near disruption and there are relatively few massive stars. As 
expected, the average mass of the stars that leave the cluster (indicated by the white line in 
Figure 6) slightly decreases throughout the inspiral (see also Fig. 7 of Giirkan et al. 2004). As 
the cluster evolves, massive stars sink towards the cluster center, leaving behind less massive 
stars which are preferentially removed. Our choice of initial conditions (requirement of core 
collapse before the end of inspiral) allows massive stars to segregate faster than the shrinking 
of the Jacobi radius because of the inspiral, so most of them remain bound throughout. 

We present the surface density of initial masses of the stars that leave the cluster and 
compare this with the assumed mass density profile of the Galaxy in Figure 8. Throughout 
the inspiral, the contribution to the background surface density is small. The only detectable 
signature of a cluster inspiral will be the stars left behind that are significantly more massive 
and hence brighter than the stars of the Galactic background distribution. 

4.3. The Cusp Retained at Disruption and IMBH 

Near disruption, the cluster forms a power-law cusp, p oc r~", around the central point 
mass. We estimated the power-law exponent to be 1.35 < a < 1.60, which is compatible 
with the results of theoretical calculations (Bahcall & Wolf 1977) and A^-body simulations 
(Preto et al. 2004; Baumgardt et al. 2004a,b) which yield a value ~ 1.55. 

The mass function of the stars in this cusp is quite different from the Salpeter IMF 
which we adopted at t = 0. The heavy part of the mass function is more populated because 
of mass segregation, even though the massive stars are preferentially removed via mergers 
with the central point mass. We show a comparison of the mass function for the innermost 
2000 stars with the Salpeter IMF in Figure 9. 

At disruption, the density of stars from the cluster exceeds the density of background 
stars. The (3-d) velocity dispersion at disruption is about 90 kms~^. Since this is much 



- 16 - 




Fig. 5. — Evolution of Model 5c (see Table 1 for initial parameters). The top panel shows 
0.01%, 50% and 95% Lagrange radii. The second panel is the growth of the central mass, 
starting at the core collapse. The third panel shows the evolution of the bound mass of the 
cluster and the bottom panel is the distance of the cluster from the GC. 




Fig. 6. — Mass and distance distribution of stars that leave the cluster during inspiral. The 
grey scale for the number of stars is logarithmic except for the ^ 1 bin; the average mass in 
each distance bin is shown by a white line. 
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Fig. 7. — Number of massive stars that left the cluster at different distances from the GC. 
The solid and dashed lines indicate the stars whose initial mass is larger than IOMq and 
6OM0, respectively. 
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Fig. 8. — Surface density of stars that leave the cluster (solid line) compared with the 
Galactic mass density used in calculations (dashed line). Note that for the cluster, the 
surface density is calculated using the initial masses of the stars. 
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Fig. 9. — Comparison of the mass function for the innermost 2000 stars at cluster disruption 
(sohd hne) with the Salpeter IMF (dashed line; normalized to the same number of stars). 
The heavy part of the mass function at disruption is much more populated because of mass 
segregation. 



- 21 - 



smaller than the velocities near the GC, the stars would stay together and be observed as 
a clump, similar to the Central cluster (Figer 2004), or the rotating disk-like structures 
observed by Genzel et al. (2003) and Horrobin et al. (2004). Additionally, the stars which 
lie very close to the IMBH will remain bound to it longer than the stars lying outside, since 
these stars will be on Keplerian orbits and their relaxation time will be longer, preventing 
the expansion of the cluster. Such a collection of stars may exhibit a structure similar to 
that of the IRS 13E cluster observed by Maillard et al. (2004). 

It is not possible to follow reliably the evolution of the cluster near disruption with our 
method, for several reasons. When the number of stars decreases to a very small value {N ~ 
1000, see Henon 1973, §2.4.1), the relaxation time becomes comparable to the dynamical 
time, making it necessary to take the effect of large-angle scatterings into account. This 
would be possible in principle, by incorporating these interactions into the Monte Carlo 
scheme in a way similar to collisions (Freitag & Benz 2002), if the rate of these events 
were low and the relaxation was still the dominant process for the evolution of the cluster. 
However, in this regime, the dynamical evolution, in particular the evaporation, of the cluster 
will not progress on a relaxation timescale and will depend sensitively on these large-angle 
scatterings. In addition, as the stars leave the cluster, the potential will get shallower. As a 
result, the restoring force on the central black hole, which keeps it near the cluster center, 
will decrease and its wandering will increase (which we do not model). 



5. Discussion 

The central parsec of our Galaxy hosts many young stars whose presence there poses 
a problem because of the unfavorable conditions for their formation. The range of age 
estimates and the peculiar dynamical properties (disk-like structures, clustering, possibility 
of larger than normal eccentricities) suggest that there is more than one mechanism at play. 
In this paper we consider one of the possibilities. 

We investigated whether it is possible to explain the presence of the young massive stars 
near the GC by bringing them there as members of a cluster. At the same time, we want 
this cluster to form an IMBH, which can then make the stars in its vicinity migrate to inner 
orbits very near Sgr A*. This scenario requires: (1) rapid inspiral of the cluster, tin ^ 3 Myr; 
(2) the disruption to take place close to the GC, R^is ^ 1 pc; (3) the cluster to undergo core 
collapse during the inspiral to form an IMBH. Our work extends that of Gerhard (2001) and 
McMillan & Portegies Zwart (2003), who carried out semi-analytical calculations, and of 
Portegies Zwart et al. (2003), Kim & Morris (2003) and Kim et al. (2004), who used iV-body 
simulations of dense clusters spiraling into the GC. 
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To determine the initial conditions suitable for this scenario, we carried out semi- 
analytical calculations and dynamical Monte Carlo simulations. A comparison of the results 
obtained by these methods showed that the change in the cluster structure has a considerable 
influence on the disruption distance. Following the evolution of the cluster also allowed us 
to draw conclusions about the demographics of the stars left behind and brought into the 
central parsec. The expansion of the outer parts of the cluster because of the relaxation 
enhances the mass loss and leads to earlier disruption. 

Our simulations show that for clusters with Rq > 10 pc, an initial mass of at least 
10^ Mq is necessary. This is in agreement with the findings of Kim & Morris (2003). We 
also found that high concentrations (Wq > 8) are required to undergo core collapse during 
the inspiral. However, if the core collapse requirement is relaxed, clusters with moderate 
concentrations {Wq ~ 6) can also survive the inspiral down to i?dis < 1 pc, starting with 
Rq ^ 10 pc and niQ > 10^ M©. This is because core collapse requires significant relaxation 
and clusters with low concentration cannot survive the accompanying expansion. 

The density of the stars that leave the cluster during the inspiral is generally small 
compared to the Galactic background density, except near disruption. Almost all the massive 
stars that leave the cluster during the inspiral do so near disruption. As a result of the mass 
segregation, most of the cluster mass close to the IMBH is in heavy (m > 10 Mq) stars. Upon 
disruption, these stars will end up on orbits close to the IMBH. As a result of their proximity, 
they can then undergo strong interactions with the IMBH and possibly get scattered into 
orbits closer to the GC. If such clusters harboring IMBHs regularly form in the GC region, 
the central parsec may be hosting more than one IMBH and this process can be realized by 
participation of multiple IMBHs. The age spread of the young stars in the central parsec 
(~3-10Myr) implies that there has been more than one instance of star formation in the 
recent history of the GC region. 

Although the discussion in this paper has focused on our Galactic center, these ideas 
also have important consequences for other galaxies and for extragalactic astrophysics. The 
direct injection into the center of a galaxy of many IMBHs produced by coUisional runaways 
in nearby young star clusters provides a potential new channel for building up the mass of a 
central SBH through massive BH mergers (Portegies Zwart & McMillan 2002). In contrast, 
minor mergers of galaxies are unlikely to produce massive BH mergers, as the smaller BH 
will rarely experience enough dynamical friction to spiral in all the way to the center of 
the more massive galaxy (Volonteri, Haardt, & Madau 2003). Our scenario has important 
consequences for LISA, the Laser-Interferometer Space Antenna, since the inspiral of an 
IMBH into a SBH provides an opportunity for direct study of strong field gravity and for 
testing general relativity (Collins & Hughes 2004; Culter & Thorne 2002; Phinney 2003). 
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Although the SBHs found in bright quasars and many nearby galactic nuclei are thought to 
have grown mainly by gas accretion (e.g. Fabian & Iwasawa 1999; Haehnelt, Natarajan, & 
Rees 1998; Richstone 2004; Soltan 1982), current models suggest that LISA will probe most 
efficiently a cosmological massive BH population of lower mass, which is largely undetected 
(Menou 2003). LISA will measure their masses with exquisite accuracy, and their mass 
spectrum will constrain formation scenarios for high-redshift, low-mass galaxies and, more 
generally, hierarchical models of galaxy formation (e.g. Haehnelt & Kauffmann 2002; Hughes 
& Holz 2003; Sesana, Haardt, Madau, & Volonteri 2004; Volonteri, Haardt, & Madau 2003). 

The main difficulty with our scenario is the requirement of large initial cluster masses 
for rapid inspiral from Rq > 10 pc. A massive (mo > 10^ Mq) cluster with a Salpeter IMF 
will contain a larger number of massive stars than currently observed at the GC region, both 
within and outside the central parsec. In particular, outside the central parsec, only a few 
young stars have been observed (Cotera et al. 1999), whereas we predict a larger number of 
these stars would have left the cluster during inspiral. There are multiple ways to resolve this 
problem, but it is not possible to determine by observation which, if any, of these mechanisms 
is most important. One of them is the formation of a cluster with a steeper IMF at the 
higher mass end, hence suppressing the number of massive stars. The upper mass cutoff we 
have chosen (mmax = 120 Mq) is equivalent to introducing a steep IMF (Weidner & Kroupa 
2004), but even more conservative choices are plausible (P. Kroupa, private communication). 
The requirement of a large initial mass is somewhat relaxed if the cluster is initially on an 
eccentric orbit (Kim et al. 2004). In this paper, we only considered circular orbits, so the 
mass requirements we find can be seen as upper limits. The problem of many young stars 
leaving the cluster during inspiral would be largely avoided by the presence of initial mass 
segregation, which is supported both by observations and theoretical arguments (Murray & 
Lin 1996; Bonnell & Davies 1998; Raboud & Mermilliod 1998; de Grijs et al. 2003). If the 
most massive members of the cluster start their life closer to the center, they will have less 
of a chance to leave the cluster outside the central parsec. We also note that the lack of 
observed young stars outside the central parsec could also result from higher extinction in 
this region (F. Yusef-Zadeh, private communication). 

The initial mass function of the Arches cluster is thought to be significantly shallower 
than the Salpeter mass function (Figer et al. 1999; Kim et al. 2000). If a cluster is formed 
with a such an IMF, the central mass will grow faster and larger, since there are more massive 
stars. This may decrease the required total mass of the cluster, at the expense of increasing 
the fraction of mass in massive stars, since a larger central object can support a cluster more 
efficiently. We carried out a few simulations with shallower IMFs, power-law mass functions 
with a = 2.00andl.75. To make a fair comparison, we kept the total mass and the total 
number of stars in these cluster the same as in our Salpeter IMF models, and hence modified 
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the minimum mass in the cluster. This way, the relaxation time is kept constant, and the 
rate of expansion is expected to be similar for different IMFs. We found that, even though 
central masses up to twice as large were grown, the changes in the disruption distance were 
very small in these simulations. We conclude that growing a large enough central mass to 
have significant impact on this scenario ( 10% of initial total cluster mass Kim et al. 2004) 
is not possible by a modification of the IMF. Note that the mass function of the innermost 
2000 stars at cluster disruption (Fig. 9), a result of the mass segregation, is very similar 
to what is observed in the Arches cluster (Stolte et al. 2002, Fig. 14). This indicates that 
the shallowness of the current mass function of the Arches can be explained by dynamical 
evolution, without the need for a shallow IMF. 

The young compact cluster IRS 13E observed near the GC region provides strong sup- 
port for the cluster inspiral scenario. This cluster has a projected diameter of ~ 0.5" and 
is at a projected distance of 3.6" from Sgr A*. Its members have a common proper motion 
with velocity 280kms~^ and no apparent common radial velocity, putting the cluster on 
an eccentric orbit. From these data Maillard et al. (2004) conclude that there is an IMBH 
at the center of this cluster with mass ~75O-13OOM0. It is possible that these stars are 
the central part of a larger cluster as in the five bright stars in the Quintuplet cluster, but 
iV-body simulations show that such a cluster, consisting of about 2000 stars, will rapidly 
experience core collapse {tec ~ 10^ yr; H. Baumgardt, private communication). Even if there 
is no IMBH in the IRS 13E cluster, its presence is a strong indication that clustered star 
formation took place in the last few million years in the GC region, close enough to Sgr A* 
that part of the cluster is now within the central parsec of the Galaxy. 

Alexander & Livio (2004) proposed that the young population observed very close (< 
0.05 pc) to Sgr A* are stars which reached this region on radial orbits and displaced stellar 
mass black holes that reside there. Gould & Quillen (2003) suggested that a member of 
this population, SO-2, is the remnant of a massive binary on an eccentric orbit which got 
disrupted by the supermassive black hole. These scenarios naturally complement the one we 
investigate in this paper. The stars scattered by the IMBH will end up on radial orbits and 
can further interact with and displace the stellar-mass black holes, or, if they are binaries, 
they can get disrupted while on these orbits. This will leave them on less eccentric orbits 
with smaller semimajor axes. The orbits of stars in the vicinity of an inspiraling IMBH have 
recently been investigated numerically by Levin et al. (2005), who showed that eccentric 
orbits with very small apocenter distances are possible if the inspiraling IMBH is on an 
eccentric orbit. 

Despite these numerical investigations by Levin et al. (2005), the final fate of a cluster 
with an IMBH, which disrupts within the central parsec, remains highly uncertain. More 
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theoretical work is needed to understand both the final phases of the disruption and the 
interaction of the IMBH with the surrounding stars. Further observations of the GC, in 
particular proper motion measurements of the closest stars to understand their eccentricity 
distribution, and higher resolution observations of the IRS 13E cluster to resolve its structure 
and dynamics, will put constraints on the various scenarios proposed. 
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near the GC, Holger Baumgardt for providing his simulation results prior to publication. 
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Kalemci for help with the preparation of some of the figures, and Ruadhan O 'Flanagan for 
comments on the manuscript. This work was supported by NASA Grants NAG5-13236, 
NNG04G176G and NSF Grant AST-0206276, and was finalized while MAG was at the 
Kavli Institute for Theoretical Physics at the University of California, Santa Barbara, as a 
Graduate Fellow. 



A. Derivation of Inspiral Rate Formulae 

In this appendix we derive the formulae for the inspiral rate dR/dt of a point mass, 
given a broken power- law Galactic mass distribution M{R) oc and an additional central 
point mass Msbh- We use Rb = 0.38 pc to denote the breaking radius, and the subscripts 1 
and 2 for regions R ^ R^ and R> R^ respectively. We follow and generalize the method of 
McMillan & Portegies Zwart (2003, see also Binney & Tremaine 1987 §7.1). 

For R ^ Rb, the enclosed mass for a circular orbit is given by 

M(i?) = Msbh + Aii?"^ (Al) 
The potential corresponding to this mass distribution is 

R tti — 1 

and the circular velocity is given by 

v.iRf = + G A.R'^^-K (A3) 

R 

Combining these, we obtain the energy per unit mass for a circular orbit, 

ZK Z tti — i 
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The acceleration resulting from the dynamical friction on an object of mass m is given by 
(Binney & Tremaine 1987, §7.1) 



a 



■f 



2 ____ 



-47rlnAxG pm— 



(A5) 



(see McMillan & Portegies Zwart 2003, for a description and numerical values of the con- 
stants). Hence the work done by dynamical friction per unit mass is 

.2 pm 



Wf = -4n In A Y 

^ Vc 

The time derivative of energy per unit mass in Equation (A4) is 



dE, _ fGMsBR , «i + l 
dt 



(A6) 



(A7) 



2i?2 2 

By setting this equal to work done by dynamical friction, we obtain the following expression 
for orbital decay 

dR _ -xlnAGi/Miaimi?(2"i-i)/2 
^ ~ (MsBH + + l)/2)(MsBH + Aii?-i)i/2 

or in dimensionless form 



(A8) 



d^ 



27T aix In A m Mi 



^{2ai-l)/2 



dr (MsBH + Mi)i/2 (MsBH + Mi^^lai + 1)/2)(Msbh + Mi^-^y/^ ' 
where we have used 



^_ R _ t ^ 
^ = ^5 ^ = 7^, To 



27r 



JTLq 



.1/2 



Ml = Ai/2^ 



(A9a) 



(A9b) 



G(MsBH + Mi; 

with Rq being the initial distance from the GC. For R > Rf, the equation for the enclosed 
mass takes the form 

M{R) = MsBH + Aii?r - ^2RT + ^2i?°' . (AlO) 

Hence, the above formulae can be used for this regime by simply substituting 2 for 1 in 
subscripts and replacing Msbh by (Msbh + — ^2-^^^)- However, the above derivations 

are valid for a 1. For a2 = I, which is the value indicated by observations (Genzel et al. 
2003), a slightly different expression is obtained 

d^ _ 4:71 X In A Mam ^^/^ 

d^ ~ ~ (Msbh + ^li^^ - ^2Rb + M2)i/2 (Msbh + Aiii^ " ^2i?6 + MaO^'/' 
where 



, (Alia) 
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i?0 



[G (Msbh + ^iK' - ^2i?b + Ms) J 



1/2 



M2 = Asi^o 



(Allb) 



For R given in pc, Ai = 2.46 x 10^ Mq, and ^2 = 2.18 x 10^ Mq according to values quoted 
by Genzel et al. (2003). 
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